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ABSTRACT 

Microlensing can provide an important tool to study binaries, especially those composed of faint or dark ob- 
jects. However, accurate analysis of binary-lens light curves is often hampered by the well-known degeneracy 
between close (i < 1) and wide (i > 1) binaries, which can be very severe due to an intrinsic symmetry in the 
lens equation. Here s is the normalized projected binary separation. In this paper, we propose a method that can 
resolve the close/wide degeneracy using the effect of a lens orbital motion on lensing light curves. The method 
is based on the fact that the orbital effect tends to be important for close binaries while it is negligible for wide 
binaries. We demonstrate the usefulness of the method by applying it to an actually observed binary-lens event 
MOA-2011-BLG-040/OGLE-2011-BLG-0001, which suffers from severe close/wide degeneracy. From this, 
we are able to uniquely specify that the lens is composed of K and M-type dwarfs located at ~'3.5 kpc from the 
Earth. 

Subject headings: gravitational lensing: micro - binaries: general 



1. INTRODUCTION 

Currently, about 2000 microlensing events are being de- 
tected each year from survey experiments (OGLE: Udal- 
ski 2003; MOA: Bond et al. 2001; Sumi et al. 2003). 
Among them, a considerable fraction are produced by ob- 
jects composed of two masses. When a lensing event is pro- 
duced by a binary object, the resulting light curve deviates 
from the smooth and sy mmetric form of a single-lens event 
dMao & Paczvhskill991h . Analyzing a binary-lens light curve 
is important because it enables to extract information about 
the lens. For general binary-lens events, this information in- 
cludes the binary mass ratio, q, and the projected separation 
between the lens components normalized by the angular Ein- 
stein radius 6'e of the lens, s (normalized projected separa- 
tion). For events with resolved caustic crossings, it is possible 



to measure or constrain the mass and distance to the lens. Not 
being dependent on the lens luminosity, microlensing pro- 
vides an important tool especially for the study of binaries 
composed of faint or dark objects. 

However, accurate analysis of a binary-lens light curve is 
often hampered by the ambiguity in fits of observed light 
curves. Many cases of degenerate solutions are accidental, 
implying that the degeneracy can be resolved by precise and 
dense coverage of light curves. However, due to a symmetry 
in the lens equation itself, degeneracies between close {s < 1) 
and wide is > 1) binaries can be very severe. This was first 
pointed out for a special case o f low-mass (planetary) com- 
panions by iGriest & Safazadeh* (1998). It was then gener- 
alized to all binaries by Dominik (1999) and analyzed still 
more thoroughly by lAnI (120051) . Therefore, devising a method 
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resolving the close/wide degeneracy is important to uniquely 
characterize binary lens systems. 

In this paper, we propose a method that can resolve the 
close/wide degeneracy. The method is based on the fact that 
orbital motion tends to be important for close binary lenses 
while it is negligible for wide binary lenses. We demonstrate 
the usefulness of the method in the interpretation of an ac- 
tually observed lensing event MOA-201 1-BLG-040/OGLE- 
2011-BLG-OOOl. 

2. DEGENERACY AND RESOLUTION 

The existence of an additional lens component can make 
lensing light curves greatly different from that of a single- 
lens event. The most dramatic feature of binary lensing is 
caustics, which represent positions on the source plane where 
the lensing magnification of a point source becomes infi- 
nite ( Schneider & Weifil ll986). As a result, binary-lens light 
curves resulting from caustic-crossing source trajectories ex- 
hibit characteristic spikes. Binary lensing caustics form a sin- 
gle or multiple sets of closed curves each of which is com- 
posed of concave curves that meet at cusps. The caustic topol- 
ogy is broadly classified into three categories depending on 
the binary separation (Erdl & Schneider 1993). For a resonant 
binary with a normalized separation similar to the Einstein ra- 
dius of the lens (i ^ 1), a single large caustic with six cusps 
forms around the center of mass of the binary. For a wide 
binary, there exist two four-cusp caustics each of which is lo- 
cated close to each lens components. For a close binary, there 
are three caustics with one four-cusp caustic located near the 
barycenter of the binary while the other two three-cusp caus- 
tics are located away from the center of mass. 

The close/wide degeneracy in the interpretation of a binary- 
lens light curve occurs due to the similarity in shape between 
the central c austic of a c lose binary and either caustic of a 
wide binary jAlbrow et al.. 2001). The binary-lens parameters 
of the pair of degenerate close/wide binaries are related by 



and 



?w = ?c(l-?c) 



S„ = S;, ' ( 1 + g'c)(?c^ - ?c + 1 )' 



-1/2 



(1) 



(2) 



where (ic,<?c) and (iw,<?w) represent the projected separa- 
tions and the mass ratios of the degenerate sets of t he close 
and wide binaries, respectively (lAlbrow et al.l 120021) . In the 
present context, it is important to emphasize that the degen- 
eracy becomes especially severe when the normalized sepa- 
ration is either substantially smaller (s ^ 1) or larger ^ 1) 
than the Einstein radius (Shin et al. 2012). 

The existence of a binary companion additionally affects 
lensi ng light curves due to the orbital motion of the bi- 
nary ( Albrow et al.''2 000t|Penn7etai]|20Tlt[Shin et alJlIoTll: 
[Skowron et al. 2011). The effect of the orbital motion is 
twofold. First, it makes the projected binary separation vary 
in time, causing the shape and size of a caustic to vary during 
the event. Second, it also causes the projected binary axis to 
rotate with respect to the source trajectory, which is equivalent 
to the lens-source relative motion deviating from rectilinear 

The orbital effect of a binary lens can provide a useful tool 
that allows one to resolve the close/wide degeneracy. The 
physical Einstein radius of a typical Galactic bulge event is 
^1-2 AU. Then, orbital periods of wide binaries would be 
substantially longer than a year. Considering that typical du- 
rations of binary-induced deviations are of order of 10 days, 
then, the orbital effect on the Ught curve of a wide binary event 
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Fig. 1.— Light curve of MOA-2011-BLG-040/OGLE-2011-BLG-0001. 
Also drawn are the best-fit curves based on single-lens, wide-binary, and 
close-binary lens models. We note that the close and wide-binary models 
are too degenerate to be visually distinguished. 

would be negligible. By contrast, orbital periods of close bi- 
naries can be small and thus the change of the lens position 
caused by the orbital motion during the deviating part of a 
lensing event can be important for some close binary events. 
Hence, if it is found that the orbital effect is needed to de- 
scribe an observed light curve, the close binary interpretation 
would be the correct one between two the possible close/wide 
binary interpretations. On the other hand, if the orbital motion 
is constrained to be extremely small, this would not absolutely 
rule out the close-binary solution, because the projected sep- 
aration of a binary with a large semi-major axis can be small 
when the binary axis is aligned with the line of sight. How- 
ever, because such extreme projections are a priori unlikely, 
strong upper limits on orbital motion would statistically favor 
wide solutions. 

3. DEMONSTRATION 

We apply the orbit-based method of resolving the 
close/wide binary degeneracy to an actually observed 
event MOA-201 1-BLG-040/OGLE-201 1-BLG-OOOl . The 
event was discovered from the survey conducted by 
the Microlensing Observations in Astrophysics (MOA) 
group on 2011 -Mar- 12 toward the Galactic bulge at 
(a, (5)12000 = (17''54"00!02 -29°06'05".70), i.e., (l,b) = 
(0.799°, -1.657°), based on observation using the 1.8 m tele- 
scope at Mt. John observatory in New Zealand dSumi et alj 
I2011h . The event was also independently discovered in 
May 2011 by the Optical Gravitational Lensing Experi- 
ment (OGLE), immediately after commissioning of the Early 
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Fig. 2. — Geometry of the lens systems con'esponding to the close (upper 
panel) and wide-binary (lower panel) solutions. In each panel, the cuspy 
figures indicate the caustics at the times marked in the upper panel. The dots 
marked by M\ and Mo represent the binary components, where the bigger 
dot represents the heavier component. The curve with an arrow represents 
the source trajectory. Right panels show enlargements of the region around 
the caustic where the deviation occuiTed. The coordinates (^,»y) are centered 
at the center of mass of the binary and all lengths are scaled by the angular 
Einstein radius corresponding to the total mass of the binaiy lens. The ^-axis 
is set as the binary axis at the time HJD' = HJD-2450000 = 5640.0. 

Warning System (EWS) of microlensing event detections 
dUdalsk i 2003) for OGLE-IV phase. OGLE-IV has been con- 
ducted since March 2010 with the 1.3 m Warsa'w telescope 
and 32 CCD detector mosaic camera at Las Campanas Ob- 
servatory in Chile. Reductions of data ■were processed using 
photometry codes developed by the individual groups. 

Figure |l] shows the Hght curve of MOA-201 1-BLG- 
040/OGLE-201 1-BLG-OOOl. It exhibits a deviation from the 
unperturbed single-lens light curve near the peak. The devia- 
tion is composed of multiple features. One at HJD'(= HJD- 
2450000) ^ 5639, shows a typical caustic-crossing feature of 
a strong spike. Another feature is the bump at HJD' ^ 5649. 
Located between the two features is an extended negative de- 
viation region that lasted for ^6 days. 

In modeling the light curve of the event, we search for a so- 
lution of lensing parameters that best describes the observed 
light curve. For the basic description of a binary-lens light 
curve, seven lensing parameters are needed. The first three 
describe the relative lens-source motion, including the time of 
the closest of the source star approach to the center of mass of 
the binary lens, fo, the lens-source separation at that moment. 
Mo, and the time scale for the source to cross the Einstein ra- 
dius, fE (Einstein time scale). Another three parameters are re- 
lated to the binary nature of the lens, including the normalized 
projected binary separation, s, the mass ratio between the bi- 
nary components, q, and the angle between the source trajec- 
tory and the binary axis, a. The last parameter is = 6'^/6'e, 
the angular source radius 0* normalized by the angular Ein- 
stein radius. The normalized source radius is needed to 
account for the finite-source effect which affects lensing mag- 
nifications during caustic crossings or approaches. With these 
parameters, the solution of the lensing parameters is searched 
by minimizing x' in a parameter space. This is done by a 
combination of grid search and down-hill approach. The grid 
search is done in the space of the parameters (s,q,ct), which 
are related to features of lensing curves in a complex way. 
For the minimization in the down-hill approach, we use 
the Markov Chain Monte Carlo method. Considering the pos- 
sibility of the existence of degenerate close and wide binary 
solutions, the search is done in the parameter space that is 
wide enough to encompass both binary solutions. From this 



TABLE 1 
Model Parameters 



Parameter 


i^iose 


wiae 


xVdof 


7734/7955 


7856/7955 


'o (HJD') 


5642.40±0.02 


5587.30±0.27 


"0 


-0.0914±0.0008 


2.9765±0.0084 


fE (days) 


53.14±0.39 


130.59±0.44 




0.506±0.002 


3.718±0.007 




1.502±0.013 


0.142±0.002 




-0.530±0.002 


1.447±0.001 


P. (10-') 


2.15±0.02 


0.88±0.01 




-0.041 ±0.067 


-0.065±0.003 




0.127±0.007 


0.1 85 ±0.004 


dsjdt (yr-') 


-0.347±0.027 


-2.040±0.019 


doLjdt (yr"') 


1.744±0.074 


-0.273±0.008 



Note. — HJD' = HJD-2450000. The lensing parameters fo and hq are 
measured with respect to the barycenter of the binary lens. The parameters 
.So and denote the binary separation and the source trajectory angle at 
HJD' = 5640.0, respectively. 



search, we find a pair of degenerate close and wide binary 
solutions. 

Although the overall shape of the observed light curve can 
be fitted with the basic binary lensing parameters, there are 
significant residuals to the fit. This gives a hint that second- 
order effects should be taken into account. Considering a long 
duration of the event, which lasted ^100 days, we examine 
the effect of the lens orbital motion and the parallax effect. 
The latter effect is caused by the positional change of the ob- 
server because of the Earth's orbital motion around the Sun 
(Gould 1992). It is known that the parallax effect results in a 
long-term deviation in lensing light c urves similar to th e ef- 
fect of the binary-lens orbital motion (iBatista et al.ll201 iT) and 
thus we consider both effects simultaneously. To first order 
approximation, the lens orbital motion is described by two pa- 
rameters dsjdt and da/dt, which represent the change rates 
of the normalized binary separation and the source trajectory 
angle, respectively. The parallax effect requires another two 
parameters tt^^m and tte^e, which are the components of the 
lens parallax vector tte projected on the sky along the north 
and east equatorial coordinates, respectively. The direction of 
the parallax vector corresponds to the relative lens-source mo- 
tion in the frame of the Earth at a specific time of the event. 
The amplitude of the parallax vector is the ratio of the size 
of the Earth's orbit to the Einstein radius projected on the ob- 
server plane. We find that taking account of these second- 
order effects significantly improves the fit for both the close 
and wide binary solutions. 

In Figure [T] we present the model light curves on top of 
the observed data for both the close and wide binary solu- 
tions. The lensing parameters of the corresponding solutions 
are listed in Table [T] It is found that the close binary solu- 
tion provides a better fit than the wide binary solution with 
Ax~ ^122. However, the degeneracy is quite severe and thus 
it is difficult to visually distinguish the two model curves. To 
be noted is that the obtained lensing parameters s and q of the 
pair of degenerate solutions well follow the relations in equa- 
tions (1) and (2), suggesting that the degeneracy is rooted in 
the symmetry of the lens-mapping equation. In Figure |2] we 
also present the geometry of the lens system corresponding 
to the individual solutions. For both solutions, the caustic- 
crossing feature of the light curve occurred when the source 
crossed a tip of a caustic and the bump was produced when the 
source approached a cusp of the caustic. The extended nega- 
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Fig. 3. — Distributions of Ax' in the space of lensing parameters. For better view, the distribution of the lens-parallax parameters, tze.n and ttee, are presented 
in a separate panel on the upper right comer 



tive deviation region occurred when the source passed along 
the fold of the caustic. Figure [3] shows Ax^ distributions in 
the lensing-parameter spaces presented to show uncertainties 
and correlations between lensing parameters. The presented 
distributions are for the close-binary solution, which provides 
a better fit. 

For both solutions, the orbital effect is important as indi- 
cated by the large values of the orbital parameters ds/dt and 
da/dt. According to the definition, ds/dt = LO implies that 
the binary separation changes by = 1 .0 per year Similarly, 
da/dt = 1.0 implies that the binary axis rotates by 1 radian per 
year The obtained orbital parameters are ds/dt = -0.35 and 
da/dt = L74 for the close binary model and ds/dt = -2.04 
and da/dt = -0.27 for the wide binary model. Consider- 
ing that the binary-induced deviations lasted ^15 days, the 



changes of the binary separation and the rotation angle over 
the course of lensing magnification are very significant for 
both solutions. From comparison of the model light curves 
obtained with and without considering the orbital effect, we 
find that the orbital effect is especially important for the de- 
scription of the extended negative perturbation region. 

The severity of the orbital effect strongly suggests that the 
event was produced by a close binary rather than a wide bi- 
nary. For a quantitative proof of this, we evaluate the ratio of 
transverse kinetic to potential energy 



KE\ (rx/AU)' 
PE y ^ 87r2(M/Mo) 



1 ds 
s dt 



da\" 
dt J 



(3) 



which must obey (KE/PE)_l < KE/PE, where KE/PE is 
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Fig. 4. — OGLE color-magnitude diagram of stars in the field where the 
binary-lens event MOA-2011-BLG-040/OGLE-2011-BLG-0001 occurred. 
The red and blue dots represent the centroid of the giant clump and the loca- 
tion of lensed star, respectively 

the ratio of (3-d imensional) kinetic to potential energy 
dPong et alj|2009l) . Here M is the total mass of the binary 
system, r± = sDi^Oe is the physical projected separation be- 
tween the lens components, and Dl represents the distance 
to the lens. An energy ratio greater than unity implies that a 
binary system is dynamically unbound. 

To determine the energy ratio, it is required to measure the 
total mass and distance to the lens. We determine these quan- 
tities from the measured lens parallax and the angular Einstein 
radius by 



M = 



and 



KITE 



AU 



7rE6'E + 7rs ' 



(4) 



(5) 



where tte = (TTE,N^ + '^E,E^y^^, tts = AU/Ds, andDs represents 
the distance to the source star The Einstein radius is mea- 
sured by 9e = Oi,/ Pi,. The normalized source radius p^, is de- 
termined from light curve fitting. In this process, we consider 
the limb-darkening variation of the source star's surface by 
modeling the surface-brightness profile as S\ = (F\/T:9i,^)[l — 
T\{l-3cosip/2)], where Fx is the linear limb-darkening co- 
efficient. Fa is the source star flux, and ip is the angle between 
the normal to the source star's surface and the line of sight 
toward the star We adopt = 0.53 and F/ = 0.44 based on 
the source type determined by the procedure described be- 
low. We estimate the angular source radius from the dered- 
dened color (V - I)o and brightness /n u sing an instrumental 
color-magnitude diagram (lYoo et al.ll2004 ). For the event, we 
estimate the instrumental color of the source based on linear 
regression method by using multi-band {V and /) OGLE data. 
Then, we measure the offset A[{V-I),I] = (-0.24,+1.41) of 
the source star from the centroid of bulge giant clump, whose 
dereddened position is known ind ependently, [(V-/),/]oc = 
(1.06. 14.41) (.Bensbv et all l20TTt iNataf et al.. .2012.) . This 



TABLE 2 
Physical lens parameters 



Parameter 



Close 



Wide 



M{Mq) 


4V^'-^^ 
"■^^-0.10 


1.80±0.05 


Mi (Mq) 


1.58±0.04 


Mz (Mq) 




0.22±0.01 


Oe (mas) 


1.17±0.01 


2.88±0.02 


jj, (mas yr"') 


8.04±0.02 


8.05±0.06 


Dl (kpc) 


T 54+0.06 
•'"'^-0.45 


1.45±0.03 


(KE/PE)x 


0.38t!!i 


9.84±0.58 



Note. — M: total mass of the close binary, Mj and M2: masses of 
the binary components, 9e: angular Einstein radius, fi: relative lens-source 
proper motion, D[^: distance to the lens, (KE/PE)x: transverse kinetic to 
potential energy ratio. 

yields the dereddened color and magnitude of the source star 
(V -1, 1)0.5 = (0.82, 15.79), indicating that the source is a low- 
luminosity bulge giant. Figure |4] shows the locations of the 
source star and the centroid of giant clump in the color- 
magnitude diagram constructed based on OGLE data. Once 
the dereddened V- / color of the source star is measured, it is 
translated into V-K color by u sing the V -I versus V - re- 
lations of lBessell & BrettI (Il988i) and then the angular source 
radius is estimated by using the r elation between V-ZT colors 
and angular stellar radii given bv lKervella etaTI (120041) . 

In Table |2] we list the projected kinetic to potential energy 
ratios for the close and wide binary solutions along with the 
determined physical quantities. For the close binary solution, 
this ratio [(KE/PE)j^ = 0.4] is smaller than unity, implying 
that the binary lens system is dynamically stable. By contrast, 
for the wide binary solution, this ratio [(KE/PE)j^ = 9.8] is 
much larger than unity, implying that the binary system is dy- 
namically unstable and thus the wide binary cannot be the 
solution of the observed light curve. This demonstrates that 
detections of orbital effect can be used to resolve severe bi- 
nary degeneracies. By resolving the degeneracy, the physical 
parameters of the lens are uniquely determined. The mea- 
sured masses of the binary components are 0.65 Mq and 0.43 
Mq, which correspond to K and M-type dwarfs, respectively. 
The distance to the lens is ^3.5 kpc from the Earth. 

4. SUMMARY AND CONCLUSION 

We proposed a method of resolving the close/wide degen- 
eracy in microlensing light curves using the effect of the lens 
orbital motion on light curves. The method is based on the 
fact that the orbital effect tends to be important for close bi- 
naries while it is negligible for wide binaries. We demon- 
strated the usefulness of the method by applying it to an ac- 
tually observed lensing event MOA-2011-BLG-040/OGLE- 
2011-BLG-OOOl. From this, we were able to uniquely iden- 
tify that the lens is composed of two K and M-type dwarfs 
located at ^3.5 kpc from the Earth. Considering that the or- 
bital effect becomes more important with the decrease of the 
binary separation, the proposed method would be important 
especially in resolving binaries with very close and wide sep- 
arations for which the degeneracy is very severe. 
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